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A host of dark energy models and non-standard cosmologies predict an enhanced Hubble rate in 
the early Universe: perfectly viable models, which satisfy Big Bang Nucleosynthesis (BBN), cosmic 
microwave background and general relativity tests, may nevertheless lead to enhancements of the 
Hubble rate up to many orders of magnitude. In this paper we show that strong bounds on the pre- 
BBN evolution of the Universe may be derived, under the assumption that dark matter is a thermal 
relic, by combining the dark matter relic density bound with constraints coming from the production 
of cosmic-ray antiprotons by dark matter annihilation in the Galaxy. The limits we derive can be 
sizable and apply to the Hubble rate around the temperature of dark matter decoupling. For dark 
matter masses lighter than 100 GeV, the bound on the Hubble-rate enhancement ranges from a 
factor of a few to a factor of 30, depending on the actual cosmological model, while for a mass of 
500 GeV the bound falls in the range 50-500. Uncertainties in the derivation of the bounds and 
situations where the bounds become looser are discussed. We finally discuss how these limits apply 
to some specific realizations of non-standard cosmologies: a scalar-tensor gravity model, kination 
models and a Randall-Sundrum D-brane model. 

PACS numbers: 95.35.+d,95.36.+x,98.80.-k,04.50.+h,96.50.S-,98.70.Sa,98.80.Cq 



I. INTRODUCTION 

Current cosmological and astrophysical observations 
clearly show that our Universe is dominated by two un- 
known and exotic components, dark matter and dark 
energy, whose energy densities are measured to fall in 
the following ranges (at 2a C.L.) 

0.092 < ficDM^ 2 < 0.124 (1) 



for the cold dark matter (CDM) component, responsible 
of structure formations and galactic and extragalactic dy- 
namics, and: 

0.30 < QuEh 2 < 0.46 (2) 

for the unclustered dark energy component, which is re- 
sponsible of the current accelerated expansion of the Uni- 
verse. (As usual f2 denotes the ratio between the mean 
density and the critical density and h is the Hubble con- 
stant in units of 100 km s _1 Mpc -1 .) 
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The nature of both components is unknown. A com- 
mon and appealing possibility is that dark matter is com- 
posed by elementary particles which decoupled from the 
thermal plasma in the early Universe. The dark energy 
component poses more serious problems: a possibility is 
that it is due to the presence of a scalar field, whose 
cosmological dynamics allows it to become the dominant 
component of the Universe just recently in the evoluti- 
onary history of the Universe. Most of the dark energy 
models predict that the expansion rate of the Universe 
may have been different from the one predicted by the 
standard Friedman-Robertson- Walker (FRW) model at 
very early stages. Not only dark energy models, but also 
other cosmological models, can predict an enhancement 
of the Hubble rate Q. If this occurs around the time 
when the dark matter particles decouple from the ther- 
mal bath, the change in the Hubble rate may leave its 
imprint on the relic abundance of the dark matter. This 
has been discussed in details in dark energy models based 
on scalar-tensor gravity |4| and in quintessence models 
with a kination phase 0, la, S an d for anisotropic ex- 
pansion and other models of modified expansion |], ^| . 
This effect implies that the basic CDM properties may 
depend on the specific cosmological model. Thus, infor- 
mation on the CDM particles may be used to constrain 
cosmological models and vice versa. In this paper we will 
exploit this connection between the dark matter and the 
expansion rate in order to derive bounds on cosmological 
models from observational data related to the dark mat- 
ter particles. We emphasize that the effects we are going 
to study do not arise because of a direct coupling between 
dark energy and dark matter: they are instead due to the 
effect induced by the dark energy model (or other models 
with enhanced expansion) on the decoupling of the dark 
matter particle. 

We will consider a generic Weakly Interacting Massive 
Particle (WIMP) as candidate of cold dark matter. For 
the cosmological model we consider models that lead to 
an enhancement of the Hubble expansion rate in the early 
Universe as compared to the rate in standard cosmolo- 
gy, like those in Refs. 0, |H EH- In order to be general 
in our analysis, we will consider a suitable parametriza- 
tion of the enhancement of the Hubble rate in the early 
Universe. As specific examples we will then relate our 
parametrization to scalar-tensor gravity (ST) models Q, 
to models with a kination phase Q and to a Randall- 
Sundrum D-brane model |l lj . 

The enhancement of the Hubble expansion rate in the 
early Universe can affect phenomena that are sensitive 
to the exact time at which they occur. This is the case 
for the freeze-out of a WIMP. The enhanced Hubble rate 
causes an enhanced WIMP dilution. The WIMP annihi- 
lation rate therefore cannot keep up with the expansion 



as long as in the standard case. Consequently, the WIMP 
freezes out earlier in this kind of non-standard cosmolo- 
gies, leading to a higher relic abundance. The enhance- 
ment of the WIMP relic abundance can be very dramatic, 
up to a few orders of magnitude 0, |f| . The requirement 
that the WIMP is the dominant dark matter component, 
i.e. that its relic abundance satisfies the bound of Eq. 
implies that the properties of the successful WIMP are 
dramatically changed with respect to the standard FRW 
case. Since the dark matter relic abundance, apart from 
the expansion, mainly depends on the WIMP mass and 
annihilation cross section, we can say that the cosmolo- 
gical model has left its fingerprint on the dark matter. A 
fingerprint which could turn out to be one of the most 
important signatures of these dark energy and other cos- 
mological models. 

In fact the WIMPs that in the modified scenario have 
the correct relic density possess a higher annihilation 
cross section than do WIMPs selected by standard cos- 
mology. This means that indirect detection signals are 
favoured in the modified scenario. In particular, it has 
been shown that the antiproton component in cosmic 
rays is a powcrfull tool to constrain dark matter prop- 
erties, even though it is affected by large astrophysical 
uncertainties E2- I 11 R e f- combined limits in the 
astrophysical- WIMP parameter space have been derived 
and it has been shown that antiprotons can powerfully set 
bounds on WIMPs in the tens-hundreds of GeV range. 
We are going to apply the same argument here, from 
which we will be able to derive constraints on the cos- 
mological models with enhanced Hubble expansion. We 
will show that the constraints on dark matter coming 
from antiproton indirect searches strongly constrain the 
cosmological models with enhanced Hubble rate in the 
early Universe. When one studies a specific cosmolog- 
ical model, the model can be constrained by a number 
of relevant observables, related to e.g. Big Bang Nucle- 
osynthesis (BBN) 0|, Cosmic Microwave Background 
(CMB), Large Scale structure and Supernovae EH , 
weak lensing E3 and General Relativity tests [l8j . We 
will show that the antiproton bound, under suitable con- 
dition, may result in even stronger limits. 

In Section[n]we discuss our parametrization of the en- 
hanced Hubble rate, and in Section ITTll we show its con- 
nection to specific cosmological models. In Section llVl wc 
discuss the relic density calculation in modified cosmol- 
ogy, for which some useful analytical approximations are 
given in Appendix E] Section deals with the cosmic 
antiproton signal and the bounds on the dark matter. In 
Section lVll we then derive the constraints on the cosmolo- 
gical models. Section lYlII translates our bounds to some 
specific cosmological models, namely scalar-tensor cos- 
mology, kination models and Randal -Sundrum D-brane 
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FIG. 1: The Hubble rate enhancement function A(T) as a 
function of the temperature. Notice that time is running from 
right to left. The solid red line has a slope parameter v = — 1, 
the dashed green line v = 0, the dash-dotted blue v = 1 
and the dotted purple line has v — 2. For all the curves we 
have chosen r\ = 10 2 and T rc = 10~ 2 GeV. The freeze out 
temperature of the relic WIMP has been fixed to Tf = 17.3 
GeV (a case which refers, e.g., to a particle with mass of 500 
GeV and annihilation cross section (o- ann f ) = Hn 7 GeV~ 2 ). 

cosmology. Finally, in Section IVIIII we summarize our 
main results. 

II. THE HUBBLE ENHANCEMENT FUNCTION 

Let us consider a class of cosmological models that 
posses a Hubble rate H in the early Universe enhanced 
with respect to its value in standard cosmology H st &. We 
parametrize this enhancement by means of a function 
A{T): 

H = A(T)H std at early times (3) 

H = H st d at later times (4) 

This situation occurs e.g. in scalar-tensor gravity models 
and in models with a kination phase. Also some models 
with extra dimensions lead to an enhanced Hubble rate. 
We will review some of these examples in the next Section 
to show that they can all be covered by the following 
paramctrization of the enhancement function A{T): 

^(T) = l + ,(|)\anh(^) (5) 

for temperatures T > T re and A(T) = 1 for T < T Te . 
By T re we denote the temperature at which the Hubble 



rate "re-enters" the standard rate and by TV a reference 
temperature, which we specify below. Clearly the Hubble 
parameter has to approach the General Relativity (GR) 
case at some epoch T rc > 1 McV, in order not to spoil 
the successful predictions of Big Bang Nucleosynthesis 
and the formation of the Cosmic Microwave Background 
Radiation (CMB). These two key events in the history of 
the Universe are typically among the major constraints 
on DE theories. 

In fig. [T] we show A(T) for r] = 100, T rc = 10~ 2 GeV 
and four different values of the exponent is, which we 
will use throughout the paper as reference points in our 
analysis. The figure shows that at the temperature Tf 
(which here is Tf ^S> T re ) all curves take the same value 
rj (since rj ^> 1). In our analysis, we will take Tf as the 
temperature at which the WIMP freezes out in standard 
cosmology. In a situation like the one shown in the figure, 
the parameter rj gets the meaning of enhancement factor 
of the Hubble rate at the time of WIMP freeze-out. 

The slope of the curve around Tf is mostly determined 
by the exponent v. We arbitrarily use in Eq. (J^J a hy- 
perbolic tangent to assure that all curves approach unity 
as a continuos function when the temperature of the Uni- 
verse approaches T rc . However, for analytic approxima- 
tions that may help to understand our numerical results, 
a sharp jump to the standard case at T = T rc is appro- 
priate. Changes in the slope of the re-entering phase will 
be briefly discussed in Sect. IVI CI 

III. COSMOLOGICAL MODELS WITH A 
MODIFIED EXPANSION RATE 

In this section we present a short list of interesting cos- 
mological models which lead to an enhanced expansion 
rate as compared to the standard Friedman-Robertson- 
Walker (FRW) model based on General Relativity. 

A. Scalar— tensor theories 

In a scalar-tensor theory of gravity both a metrical 
tensor g^ and a scalar field ip are involved in the de- 
scription of the gravitational interaction 0] . Because of 
this additional scalar degree of freedom, the Universe ex- 
pands differently compared to a standard FRW solution 
of the Einstein gravity equations This class of theo- 
ries can be formulated in different frames related to each 
other by a Weyl re-scaling of the metric • As pointed 
out in Ref. |2l| . a frame transformation amounts in a 
change of units and therefore the physical results cannot 
depend on the frame. A frame-independent formulation 
of the scalar-tensor cosmology is given in ■ In order to 
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deal with scalar field- independent masses and couplings, 
in the language of scalar-tensor theories this means that 
we stay in the Jordan frame. This was our approach 
in the calculation of the effects induced on the WIMP 
decoupling in scalar-tensor theories in Rcf. Q . 

The ratio between the Jordan frame expansion rate 
and the standard General Relativity expansion rate is 
given by 0: 



H 2 



A\ V ) 



id 



i - (^')76 



(6) 



where ct(<p) = d\ogA(tp)/dip, A(ip) is the Weyl factor 
relating the Jordan and Einstein frame and a prime de- 
notes derivation with respect to the logarithm of the scale 
factor. 

To find a general correspondence between our para- 
metrization in Eq. JSJ) and a scalar-tensor behaviour of 
the expansion rate is quite difficult. This is due to the 
fact that the scaling of the ratio in Eq. 10 is deeply re- 
lated with the cosmological dynamics of the scalar field. 
However, as shown in a specific example in Ref. 0, a 
numerical solution of the scaling behavior of Eq. © can 
be translated in terms of our parametrization. In the 
simplest case of a slowly varying scalar field (ip' ~ 0) 
we have H = A(tp)H a td, and in the specific example of 
we have A(ip) = A(tp(Tf))(T/T{)-°- 82 . In terms of 
our parametrization, this example has r\ = A(tp(T{)) and 
v = —0.82. Fig. 6 of Ref. 0] shows that the enhance- 
ment of the Hubble rate around the WIMP freeze-out 
may be quite sizeable, up to factors of 10 4 . We remind 
that the model studied in Ref. Q was a perfectly viable 
scalar-tensor model, since it evaded all experimental con- 
straints from BBN, CMB and gravitational probe limits, 
like the Cassini mission ^^|. The enhanced Hubble rate 
reflects in an anticipated WIMP decoupling, with an en- 
suing larger relic abundance. Ref. |j| showed that for 
such a fast re-entering of the Hubble rate on its GR be- 
haviour, it may be possible that the already-decoupled 
WIMPs start a brief phase of re-annihilation, due to the 
fact that they are still over-abundant and therefore their 
annihilation rate is larger than the GR Hubble rate. This 
phenomenon has the consequence of reducing the WIMP 
relic abundance with respect to the value it would have 
had without re-annihilation. Nevertheless, the outcome 
of this scalar-tensor model is that the WIMP current 
abundance may be larger that the GR one by up to 2- 
3 orders of magnitude. This result would then be even 
stronger if a re-annihilation phase does not occur, due to 
a slower re-entering of the Hubble rate to GR. Summa- 
rizing: perfectly viable scalar-tensor models from the 
point of view of cosmological and astrophysical observa- 
tions, may predict a strongly enhanced WIMP relic abun- 
dance for a give candidate. This was our original motiva- 



tion to try to set additional limits on these cosmologies 
by means of observables related to the dark matter sec- 
tor, namely the antiproton indirect-detection signal, as 
it will be described in the next Sections. However, our 
discussion may be applied also to other cosmological sce- 
narios with enhanced early Hubble rate, like the following 
two relevant cases. 



B. Kination 

A kination is a period in which the total energy density 
of the Universe is dominated by the kinetic term of a 
scalar field. A phase of kination is generically expected in 
quintessence models based on tracking solutions for the 
scalar field, Ref. j^. When the scalar potential V(Q) 
is negligible compared to the kinetic energy of the scalar 
field, the total energy density in the scalar field, scales 
like ~ a -6 (where a is the scale factor of the Universe). 
This means that during kination H 2 cx ptot — p$ oc a~ 6 . 
More precisely, the ratio between the expansion rate H 
during a kination period and the standard expansion rate 
H s td is given by: 



Jl 2 

H: 



Id 
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where p$/p r can be written as 



p<t> 



= m 



h(T) 



KT t ) 



9(T f ) (TV 
9(T) \T f J 



(7) 



(8) 



with 77$ = p$(T{) / p r (T{) and h e g and g e g respectively 
the entropy-density and energy-density effective degrees 
of freedom. The approximation in Eq. JS} is justified 
only in a range of temperatures where h e s and g e s do 
not change considerably with respect to their value at 
T f . 

From Eq. Q and © we see that a kination model 
can be approximated by our parametrization once the 
following values for the parameters are chosen: v = 1 
and 77 = ^/fj^. 



C. Extra Dimensions 

We refer to the RSII model ^lj . In this model a single 
3-brane is embedded in a five dimensional bulk with a 
negative five dimensional cosmological constant. Unlike 
other Extra Dimension models, here the fifth dimension 
is not compact and no orbifold boundary conditions are 
imposed. Moreover, the five dimensional metric is not 
factorizable and an exponential function of the fifth co- 
ordinate multiplies the four dimensional metric. With 
this set up, a Kaluza Klein (KK) reduction of the five 
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dimensional gravitational excitations gives rise to a spec- 
trum with a zero mode (the standard four-dimensional 
graviton) localized in the extra dimension and a conti- 
nuum of massive KK modes weakly coupled to the low 
energy states on the brane. Due to the fact that the four- 
dimensional gravitons can propagate in a confined region 
of the fifth dimension and since at low energy the massive 
KK modes are only weakly coupled, the theory is not in 
contrast with experimental gravity even if the volume of 
the extra dimension is infinite. The Einstein equation on 
the 3-brane are studied in Ref. [3 and the cosmology 
of such a model is carefully investigated in Ref. |2{| . In 
particular it has been shown that the ratio between the 
expansion rate H of such a model and the expansion rate 
H st d of standard cosmology is given by: 



H 2 



Id 



P_ 

2A 



(9) 



where A is the tension on the brane and p the energy 
density on the brane. When the total energy density p 
equals the radiation energy density p r we have: 



H 



H. 
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Pr(Tf) 

2A 



(10) 



Comparing this expression with Eq. (0), we see that such 
a model is described by our parametrization with the 
values of the parameters v = 2 and rj = y' p r (T{)/(2\). 

IV. THE RELIC DENSITY CALCULATION 

In this section we will briefly discuss how a cosmolo- 
gical model with enhanced Hubble rate affects the relic 
abundance of WIMP dark matter. This effect was stu- 
died in details in Ref. for the scalar-tensor case and in 
Refs. 0,0,1113 for the kination case. See also Refs. HE3 
for anisotropic expansion and other models of modified 
expansion. 

The evolution of the WIMP number density n as a 
function of cosmological time t is described by the stan- 
dard Boltzmann equation, the only difference being that 
the standard Hubble parameter iJ st d(T) is now replaced 
by the modified Hubble rate, H(T) = A(T)H std (T): 



du 



= -3Hn - 



"eq 



(11) 



where (<7 a nn w ) is the usual thermally averaged value of 
the WIMP annihilation cross section times the relative 
velocity. The modification of the Hubble rate can be 
rephrased as a change in the effective number of degrees 
of freedom from g cff (T) to A 2 (T)g cfi (T). The Boltzmann 
equation is more conveniently solved by rewriting it in 



terms of the comoving abundance Y = n/s where s is 
the entropy density s = (27r 2 /45) h(T)T 3 and studying 
the evolution as a function of the temperature T: 



dY 



1 



^=--T>an„«)(y 2 -^ q ) 



xH 



(12) 



where x ~ m x /T. In our analysis we solve the Boltz- 
mann equation Eq. I|12|l numerically down to the current 
value of the comoving abundance Yq. The WIMP relic 
abundance is then simply: 



tt x h 2 



m x soY 



(13) 



where sq and p c are the current values for the entropy 
density and the critical mass-density of the Universe. In 
order to get some insight in our numerical results, we 
report in Appendix ^ some useful analytical approxi- 
mations which are valid for large enhancements and re- 
entering temperature (much) lower than the temperature 
at which decoupling occurs. Notice that in our definition 
of the function A(T) in Eq. JSJ, we use as a normaliza- 
tion temperature the freeze-out temperature obtained in 
standard cosmology and defined in Appendix ^] We re- 
port it here for convenience: 



X{ = In 



0.038 m p \ g m x (a a 



V -1/2 -1/2, 



(x f ) (14) 



The effect of the enhanced Hubble rate on the WIMP 
relic density can be very important. As was found in Ref. 
0, and as we are also going to see later in this paper, 
the relic density can be up to few orders of magnitude 
larger in the modified scenario as compared to the stan- 
dard case. From combined cosmological observations we 
have the very stringent bound of Eq. Q for the cold dark 
matter density. This means that the WIMP candidates 
selected in the case of enhanced Hubble rate will be vastly 
different (i.e. have different values of their relevant pa- 
rameters, like mass and couplings) from the WIMPs that 
fit into the standard cosmology picture. The cosmologi- 
cal model has, in other words, left its signature on the 
dark matter. A signature which might turn out to give 
valuable clues. 

Not only can we say that the WIMPs selected in the 
standard and modified cases are different, we are also 
able to say something general about the difference of 
their phenomenology. This builds on the fact that the 
standard WIMP relic abundance in general is approxi- 
mately inversely proportional to the WIMP annihilation 
cross section. Analytically one finds (see Appendix 
the well known behaviour: 
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(15) 
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where (c ann u)i n t is the following integration of the ther- 
mal annihilation cross section: 



(fannU)int 



1 



Q(x) (a- ann v) 



G{xi) J Xf A(x) x 2 



dx 



(16) 



The WIMPs that fulfill the relic density constraint of 
Eq. 0} when calculated with the enhanced Hubble rate, 
would have a much lower density when recalculated in 
standard cosmology. From Eqs. (|15I16|) this then means 
that these WIMPs have a higher annihilation cross sec- 
tion than WIMPs that satisfy the relic density con- 
straints when calculated with the standard Hubble ex- 
pansion. A high WIMP annihilation cross section in the 
early Universe does in general mean that also the current 
WIMP annihilation rate in the Galaxy is high [4l|. The 
WIMP candidates selected by the models of enhanced ex- 
pansion rate are therefore in general more suitable for in- 
direct detection than are the WIMP candidates selected 
in the standard scenario. At the same time this means 
that the bounds on the WIMP annihilation cross section 
coming from searches for indirect WIMP signals can be 
used to constrain the expansion rate in the early Uni- 
verse. In this paper we are going to analyze the con- 
straints coming from the cosmic antiproton signal, which 
is the only indirect probe for which strong constraints 
can be determined 



V. THE COSMIC ANTIPROTON SIGNAL AND 
THE BOUND ON THE ANNIHILATION CROSS 
SECTION 



The antiproton component of cosmic rays has been 
measured by many detectors in space and the most recent 
results come from BESS, AMS and CAPRICE experi- 
ments. The standard production of cosmic antiprotons 
from spallation of nuclei on the diffuse Milky-Way gas 
is enough to explain the data |2(| but the error- 
bars leave a small room also for an exotic antiproton 
signal. Such a signal could come from the annihilation of 
WIMPs in the Galaxy (see e.g. Ref. [l2( )■ This situation 
is shown in Fig. [3 where the experimental data arc plot- 
ted together with the background component and some 
examples of a signal coming from WIMP annihilation, as 
calculated in Ref. flil ]. 

The antiproton flux produced in a point of cylindrical 
coordinates (r, z) in the galactic halo, where the dark 
matter density is p(r,z), depends on the annihilation 
cross section of the WIMPs averaged over their veloci- 
ty distribution in the galactic halo (cr ann z;)o and on the 
final-state branching-ratios into the different possible fi- 
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FIG. 2: Primary TO A antiproton fluxes as a function of the 
antiproton kinetic energy. The dashed line denotes the secon- 
dary component, due to spallation of cosmic rays (i.e. the 
background) taken from Refs. [irl l27|l. The other lines are 
representative fluxes from neutralino annihilation of differ- 
ent masses [l^l : the solid line refers to m x = 60 GeV, the 
long-dashed line to m x — 100 GeV, the short-dashed line to 
m x = 300 GeV and the dotted line to m x = 500 GeV. The 
astrophysical parameters for galactic propagation assume the 
best-fit values according to the analysis of Refs. [l2l l2s| . 
Solar modulation is calculated for a period of minimal solar 
activity. Full circles show the BESS 1995-97 data 0; the 
open squares show the BESS 1998 data |3oT| : the stars show 
the AMS data |3*H and the empty circles show the caprice 
data EH. 



nal sates F : 



q p {r,z : T p ) 
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P 2 (r,z) 
2m x 



^BR( X x^F) 



dNl) 
dT p J 

^ (17) 

where Tp is the antiproton kinetic energy. The antipro- 
tons then propagate and diffuse in the galactic medium 
until they reach the Earth position in the Galaxy of co- 
ordinates (Rq, 0): 



q p [r, z : T p ) — ► $p(i? Q ,0 : T p ) 



(18) 



The process of propagation has been discussed in details 
in Ref. ^3 1 where it has been shown that, contrary to the 
case of the spallation antiprotons |2s| , a large uncertainty 
in the low energy tail of the antiproton signal is present, 
due to the current uncertainties on the astrophysical pa- 
rameters which describe the diffusion and propagation 
processes 0, • This uncertainty is about one order 
of magnitude up or down around the best fit result |l2| , 
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FIG. 3: The solid line shows the maximum value (upper 
bound) of {er ann i;)o that can be allowed by the cosmic an- 
tiproton data in the fist energy-bin of BESS and AMS 
(T TOA = 0.23 GeV.) (<T ann ?;)o is the thermal average in the 
galactic halo of the WIMP annihilation cross section times 
relative velocity. The upper bound is shown as a function of 
the WIMP mass m x and has been derived using a median set 
of astrophysical parameters for the antiproton propagation. 
The upper and lower dotted lines show the uncertainty in the 
determination of the upper limit on {a ann v)o coming from the 
astrophysical uncertainties. The horizontal solid line shows 
the lower bound on {a ann v)o coming from the relic density 
of cold dark matter in standard cosmology, under the hypo- 
thesis that the WIMP annihilation cross-section in the early 
Universe (cr ann i>) equals (<T an ni , )o (like in the case of s-wave 
dominance) . 



for antiproton kinetic energies below 1 GeV where a sig- 
nal may be more promisingly approached for WIMPs in 
the tens-hundreds of GeV mass range, and without the 
need of boosted enhancements on the dark matter den- 
sity 0, . The additional solar modulation effect in- 
troduces another source of uncertainty which is expected 
to be less important than the one due to galactic diffu- 
sion. Also the uncertainty coming from the dark matter 
profile is not very relevant, since antiprotons which reach 
the Earth are produced relatively close in the Galaxy, due 
to diffusion, and therefore they do not strongly feel the 
quite uncertain galactic-center mass distribution. Dif- 
ferences in the galactic halo shape affect the antiproton 
signal at most 20% 01- For a thorough discussion of 
all these topics and of the calculation of the antiproton 
signal which we also use in the present analysis, we refer 
toRef. [13. 

The antiproton signal is a powerful tool for constrai- 



ning WIMP dark matter, since the amount of antiprotons 
produced by typical WIMPs which can account for the 
dark matter abundance of Eq. Q is at the level of the 
background and may even exceed it sizably 

E3E3- The 

astrophysical uncertainties limit somehow the capabili- 
ties of this type of signal, but nevertheless combined lim- 
its in the astrophysical-particle physics parameter space 
may be set, especially for relatively light WIMPs. Ref. 
|l3| showed how WIMPs in the mass range from few GeV 
to hundreds of GeV are currently constrained by antipro- 
tons searches. The possibility to set bounds mainly re- 
lies on the fact that the low-energy tail of the predicted 
antiproton signal may exceed the room left in the ex- 
perimental uncertainty of the data over the calculated 
background. We make use of the same argument in the 
present discussion to use the antiproton signal to set lim- 
its on (<T ann w)o and to transform these limits on bounds 
on the enhancement of the Hubble rate around WIMP 
freeze-out under the condition that the WIMPs satisfy 
the cosmological bound on the amount of cold dark mat- 
ter. 

A few comments are in order here before we proceed to 
discuss our strategy. First of all, we explain why we do 
not use also other indirect detection signals other than 
antiprotons. WIMP annihilation may obviously also pro- 
duce gamma rays and positrons. However, contrary to 
the case of antiprotons which, almost independently of 
the halo density profile, are naturally at the level of the 
background and experimental data for cosmologically rel- 
evant WIMPs, both positrons and gamma rays usually 
require sizeable boosts in the dark matter density in or- 
der to reach detectable levels. This introduces a strong 
model-dependent variable which does not allow us to use 
these indirect signals as reliable tools for setting limits. 
For instance, a limit obtained from gamma-rays would 
fade out completely unless a very steep density profile 
is present at the galactic center. The same occurs for 
positrons, which need strongly clumped structures very 
close to our position in the Galaxy. Even though these 
indirect signals are very appealing for dark matter stu- 
dies, they do not prove to be useful in the analysis we 
want to carry on in this paper. Very promising will po- 
tentially be antidcutcrons |34j . but for that we have to 
wait for the foreseen experimental set-ups able to access 
the required sensitivities ji^, |3(| . 

For similar reasons, we are concentrating our analy- 
sis on the low-energy tail of the antiproton flux: also 
at energies in the tens of GeV range (where CAPRICE 
data are available) a signal could manifest itself above 
the background, which is here fast decaying. Moreover, 
astrophysical uncertainties are less relevant at these ener- 
gies. However, in order to have large signals in this range 
of energies, we need a suitable WIMP number density, 



8 



which can be likely obtained only with some degree of 
over-density (due to the fact that we need here heavier 
WIMPs, whose number density is damped by the m~ 2 
factor) |37| . Again, since we do not want to add addi- 
tional arbitrary inputs in our analysis (the boost factor) , 
we focus on the low-energy tail. 

Finally, we must remind that the antiproton signal 
depends on (and therefore can be constrained to set 
bounds on) (a ann v)o, while the relic abundance depends 
on (<7 aml v), i.e. on the thermal average on the annihi- 
lation cross section at a different (larger) temperature. 
A constraint on (cr an ni ; )o is not therefore directly trans- 
ferable to a bound on the WIMP relic abundance, and 
vice-versa. The usual expansion of the thermal average: 

(CTannW) = (!+- + ■■■ (19) 

holds in most of the cases, noticeable differences are 
when co-annihilation effects are present [3^, |3j| or when 
annihilation occurs close to a resonance or a threshold 
[iof . However, in most of the typical cases, the ther- 
mal average of the annihilation cross section is mildly 
temperature-dependent. We will first discuss through 
the paper the case of temperature-independence, i.e.: 
(a ann v) = (c ann u)o = a. We will then discuss how our 
results are changed when the first-order expansion of Eq. 
(|19|l is relevant. The most important effect on our results 
basically comes from the difference between the freeze- 
out and current value of (<7 ann f)i *- e - on the ratio: 

K = (gaim ^ Tf (20) 

Therefore, an estimate of how our results would change, 
for instance when the relic abundance is determined by 
co-annihilation effects, is to assume a given value for 1Z 
in relating (a ann v) at freeze-out and (cr an n w )o- We also 
need to comment that, in order to apply the antiproton 
bound, the WIMP annihilation must proceed sizably to 
non-leptonic final states. In fact, for annihilation into 
leptons, antiprotons are not produced and our bounds 
are loosened by a factor given by the branching ratio 
into non-leptonic final states. 

Let us now determine the upper limit on the WIMP 
annihilation cross section derived from the observational 
upper limit on the exotic cosmic antiproton flux. We do 
not attempt here a statistical analysis of the data as was 
done in Ref. 01 . We instead use a simpler approach 
which uses the most relevant information coming from 
the antiproton data and at the same time allows us to 
have an insight on the results, which could instead be 
more difficult to obtain by a statistical treatment of the 
data. This approach was the one also adopted in Ref. 



|12(. We consider the experimental result in the lowest- 
energy bin at Tj OA = 0.23 GcV. By subtracting the 
background j2||, we are left with a 90% C.L. upper limit 
for the exotic antiproton component at the top of the at- 
mosphere (TOA) of: $J OA = 2 ■ lO-^m^s^sr^GcV- 1 
|l2j . In order to find the corresponding upper bound on 
the WIMP annihilation cross section we must specify our 
assumptions in the calculation of Eq. I|17I18|) of the ex- 
pected antiproton flux for any given mass m x and cross 
section (<r ann v) . 

We assume an isothermal halo model with core radius 
3.5 kpc and local dark matter density 0.3 GeV cnr 3 . 
When nothing else is mentioned we use the mean values 
for the set of astrophysical parameters of the propagation 
and diffusion models given in Ref.0|, but we discuss 
also the set of parameters which provide the maximal 
and minimal antiproton flux, in order to properly take 
into account this intrinsic and dominant source of un- 
certainty. As for the antiproton spectrum, we assume for 
definiteness that the WIMP annihilate only into a bb pair, 
and we therefore use the corresponding spectra given in 
Ref. A more generic annihilation final state would 

only mildly change our results (except for the already- 
mentioned leptonic final state). 

The result on the upper bound on (<7 ann i>)o is shown in 
Fig. |21 as a function of the WIMP mass and for the me- 
dian, as well as maximal and minimal set of astrophysical 
parameters. We notice that the astrophysical uncertainty 
is severe, but nevertheless allows us to set limits on the 
maximal amount of (<7 ann v)o which is allowed in order not 
to go in conflict with the experimental data. The hori- 
zontal solid line instead shows the lower limit on (cr ann f ) 
coming from the cosmological bound on the WIMP relic 
abundance of Eq. Q), in the case of temperature inde- 
pendent thermal average and standard cosmology. The 
figure clearly shows that there is tension between the cos- 
mological bound and the antiproton limit: in the case of 
the median bound from antiproton, we see that WIMPs 
lighter than about 50-60 GeV lead to an antiproton sig- 
nal which is in excess of the experimental bound. This 
tension is released if we consider the minimal set of astro- 
physical parameters (upper dotted curve) , but the figure 
clearly shows that antiprotons arc a powerful tool for 
setting limits. Ref. 0] already discussed in details this 
tensions and the corresponding limits. We instead here 
use this argument to set bounds on cosmological models: 
since the cosmological models we are discussing predict 
larger lower-bounds on (a ann v), we see that the tension 
with antiproton data is enhanced and limits can be set 
on the maximal amount of enhancement of the Hubble 
rate at the time of WIMP freeze-out. 
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VI. CONSTRAINING THE HUBBLE RATE 
WITH DARK MATTER 

In this Section wc show how the combined constraints 
on the WIMP relic density and antiproton signal can 
strongly constrain the Hubble rate in the early Universe. 
Let us first recall our assumptions about the Hubble rate 
and the WIMPs. We assume that the Hubble rate is en- 
hanced by a temperature-dependent factor A(T) in the 
early Universe, for which we consider the parametrization 
of Eq. ©. At the temperature T f at which the WIMP 
would freeze out in the standard case, the enhancement 
factor is equal to the parameter 77, A(T{) = 77 (for 77 S> 1). 
The slope of A(T) is set by the parameter v. At a tem- 
perature T = T 10 , the expansion rate "re-enters" the 
standard one. Wc consider a general WIMP model cha- 
racterized solely by a mass, m x , and an annihilation cross 
section. We assume first that (a ann v) = (a ann v) is tem- 
perature independent, and therefore it determines both 
the relic abundance and the antiproton signal directly. 
The WIMP should explain the observational amount of 
cold dark matter of Eq. ifTJl. i.e. it must be the dominant 
CDM component. 

In the following wc combine the relic density and an- 
tiproton constraints on the WIMP to obtain the upper 
bound on the enhancement of the Hubble rate. As we 
have discussed, this is possible because the Hubble rate 
together with the WIMP mass and cross section deter- 
mines the WIMP abundance, and at the same time the 
WIMP cross section is bounded from above by the cos- 
mic antiproton data. In Section IVI Al we find the mass 
dependent upper bound on 77, i.e. the enhancement fac- 
tor at the time of the WIMP freeze out. In Section fVl Bl 
we determine the corresponding upper bound on the en- 
hancement of the relic WIMP density as compared to 
standard cosmology. Finally, in Section lYl O we analyze 
what happens to the bound on 77 when some of our as- 
sumptions about the Hubble rate and the WIMP model 
are changed. 

A. Constraining the expansion rate 

Let us in this section show our numerical results on 
the Hubble rate constraints derived from the relic den- 
sity and antiproton bound. Our parametrization of the 
Hubble enhancement function Eq. (JSJ has three free pa- 
rameters: 77, v and T re (while Tf is set by the WIMP 
mass and cross section through Eq. p4[lh As we shall 
see at the end of this section, the enhancement is basi- 
cally independent of T re , as long as T re <C Tf. This can be 
understood also by means of the analytical approxima- 
tions given in Appendix lAl Wc choose to set T ro = 1CP 3 
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FIG. 4: Exclusion plot in the (77, v) plane for a WIMP mass 
of m x = 500 GeV. 77 and v are two of the parameters of the 
Hubble enhancement function. 77 is the enhancement factor at 
the freeze-out temperature of the WIMP. The "re-entering" 
temperature has been fixed at T ro = 10" 3 GeV. The WIMP 
annihilation cross section (a^nnv) is chosen for each point in 
order to fulfill the relic density constraint. This cross section 
has then been compared to the antiproton upper bound on the 
cross section. The central almost vertical band denotes points 
which have exactly the limiting value of the cross section; the 
points on the left of the band refer to cross sections below the 
antiproton upper bound; the points on the right of the band 
refer to cross section in excess of the antiproton upper bound. 

GeV throughout the paper as a reference value. This cor- 
responds to the lowest value which can be safely assumed 
to be compatible with big bang nucleosynthesis Q . The 
enhancement of the Hubble rate is then studied in the 
two-dimensional parameter space (j), v). 

Let us now construct an exclusion plot in the (77, v) 
plane. We first fix the WIMP mass at m x = 500 GeV as 
a reference value. Then, for each point in the (77, v) plane 
we make a scan in the WIMP annihilation cross section 
(a ann v). The cross section together with the mass fixes 
Tf, which is the temperature at which the WIMP would 
freeze out in the standard cosmology scenario. We then 
solve the Boltzmann equation with the enhanced Hub- 
ble rate. We can now apply the first bound, namely 
the constraint Eq. ^ on the relic dark matter density. 
This constraint selects a small interval in (cr a nni ; ) for each 
(77, v). Finally we apply the constraint coming from the 
cosmic antiproton data on (cr ann 7j)o, by means of the re- 
sult shown in FigOI which gives us the upper bound on 
the annihilation cross section for any given mass. 

In fig. 01 we show the exclusion plot in the (77, v) plane. 



The points in the central almost -vertical band refer to 
annihilation cross sections (<r ann i>) equal to the maximum 
value allowed by the antiproton data (using the median 
set of astrophysical parameters) as shown in Fig. |3J The 
points on the left of the vertical band all refer to smaller 
values of the cross section and thus satisfy the antiproton 
bound. On the contrary, the points on the right are in 
disagreement with the antiproton data as they require a 
high cross section in order to fulfill the density constraint. 
Notice that the cross section is not uniquely determined 
by the density constraint, since for fixed values of 77, v 
and m x we can slightly vary the annihilation cross sec- 
tion to move around in the allowed density interval of 
Eq. (Q. Similarly, points next to each other can have 
identical cross section and fulfill the density constraints 
with just slightly different values of the relic abundance 
[4^ . Thus, the width of the vertical band, that marks 
the antiproton bound, is due to the allowed interval of 
the relic abundance. Also, in the vertical band we do 
not only have models with a cross section equal to the 
antiproton limit, but also models with a slightly lower or 
higher cross section. 

From fig. 0] we conclude that the constraint coming 
from the cosmic antiprotons produced by dark matter 
annihilation can be used to put constraints on the Hub- 
ble rate in the very early Universe. The figure shows 
that, for m x = 500 GcV, we get an upper bound on the 
enhancement factor 77 of the Hubble rate at the WIMP 
freeze-out of around 10 2 (with some dependence on the 
actual values of the slope parameter v). This bound on 
77 may be considered as significant: in Ref. Q we found 
perfectly viable scalar-tensor models with enhancement 
factors of up to 10 4 at the WIMP freeze-out tempera- 
ture, once the BBN, CMB and gravitational probe limits 
were satisfied. Our results shows that alternative, and 
even stronger, bounds on dark energy models can there- 
fore be obtained also by looking at the imprints left on 
the dark matter properties by these modified cosmologies 
at the time dark matter formed in the early Universe, if 
we assume that dark matter is provided by a thermal 
relic. 

The results in fig. 21 were calculated for a fixed WIMP 
mass of 500 GeV. Let us now explore how the bound on 77 
depends on the WIMP mass. The result is shown in fig.|5J 
Each vertical band corresponds to the maximum allowed 
cross section for any given mass according to fig. El Going 
from left to right in fig. [SJ the vertical bands corresponds 
to the masses 60 GeV, 100 GeV, 200 GeV, 400 GeV, 600 
GeV and 1000 GeV respectively. For any given mass, 
the part of the (77, v) plane to the right of the vertical 
band is excluded by the antiproton data while the part 
to the left is allowed. The constraint on 77 is very strong, 
especially for low WIMP masses. In particular, the figure 
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FIG. 5: The same as Fig. 21 but including the result for 
different WIMP masses. Going from left to right the ver- 
tical bands correspond to the antiproton bound for WIMP 
masses of 60 GeV, 100 GeV, 200 GeV, 400 GeV, 600 GeV 
and 1000 GeV respectively. For any given mass, the verti- 
cal band marks the limit between the excluded part of the 
parameter space to the right and the part to the left which 
is in agreement with the antiproton data. The "re-entering" 
temperature is T IC = 10~ 3 GeV. 

shows that WIMPs lighter than about 60 GeV are close 
to being excluded. A word of caution is in order here: 
to make a precise claim on a lower bound on the WIMP 
mass would require a more thorough analysis, since we 
should include the uncertainty on the propagation and 
diffusion parameters for the calculation of the antiproton 
flux and also apply a more refined statistical analysis 
on the antiproton data. We will show the effect coming 
from the astrophysical parameters in Sect. IVI Gl and we 
anticipate that for the less stringent set of parameters, 
WIMPs lighter than 60 GeV are perfectly viable. These 
results are in agreement with what is found in Ref. [l^ | 
where a careful analysis of the bounds that can be set 
to the WIMP parameters by antiproton data in standard 
cosmology has been performed. 

From fig. [S] we sec that the upper bound on 77 increases 
less than an order of magnitude when the parameter v 
is varied from —1 to +1 for any given mass. For large 
masses, the increase can be two orders of magnitude when 
we continue to v = 2. This can be seen in fig. which 
displays the mass dependence of the upper bound on rj, 
as derived for four different values of v. 

The behaviour of the upper bound on 77 can be under- 
stood also by means of the analytic approximations of 
Appendix The conditions we want to satisfy at the 
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FIG. 6: Upper limit on the value of r\ obtained by imposing 
the combined constraints on the WIMP relic abundance and 
the WIMP antiproton signal. rj is the Hubble enhancement 
factor at the time where the WIMP freezes out in standard 
cosmology. The bound on rj is shown as a function of the 
WIMP mass for the slope parameter of the enhancement func- 
tion v = —1, 0, 1, 2. The parameter T rc has been set to 10 -3 
GeV. Tf is determined by the WIMP mass and cross section 
by means of the freeze-out condition of Eq. 1141 . The WIMP 
cross section is fixed at the limit of Fig. |2] 

same time are the following: 

cosmological bound (21) 
ann^)o antiproton bound (22) 

where (£lh )cdm is the maximal allowed value of the 
relic abundance constraint of Eq. Q and (<7 a nn^)o is the 
upper limit on the annihilation cross section coming from 
Fig. El This limit may be approximated as: (o'ann'y)o — 



0.95 • 10" 12 ml 95 GeV- 2 for 



bounds transform in the following condition 



> 60 GeV. The two 



fM < (<7ann«)o {^h 2 ) C DM 



(23) 



where, by means of the approximations given in Ap- 
pendix [S] 



m = c 



I(x[,x,r),is) +r G - 
x 



(24) 



and C = so Pc 1 -B G^ 1 (x/). The maximal value of r\ 
occurs when rj saturates the bound of Eq. H23|) . It is 
easy to show that: 

)o {Qh 2 ) C DM 



Vn 



C(v — 1) xf 

m v)f, (flh 2 ) C DM 



C{\ — v) X{ 



1.95 



' for v > 1 



for v < 1 (25) 
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FIG. 7: Upper bound on the Hubble enhancement parame- 
ter 77 as a function of T re , i.e. the temperature at which the 
Hubble expansion "re-enters" the standard evolution. The 
upper bound on rj is derived from the combination of the con- 
straint on the WIMP relic density and antiproton signal. The 
numerical derivation is done for two different WIMP masses 
and four different values of the parameter v. The masses are 
m x = 200 GeV and m x = 1 TeV corresponding to T t ~ 7.46 
GeV respectively Tf ~ 31.86 GeV for the limiting cross sec- 
tion. The upper set of curves is for the high mass. For each 
of the two set of curves we have v = 2, 1,0, —1 when going 
from the top to the bottom. 



and a slightly more complicated expression for v = 1 
(kination case). These behaviours are confirmed by 
our numerical calculations, where we approximately find 
?7max oc to 2 / 1 for v = —1,0,1 (with a slight depen- 
dence on which masses are used for the calculation) and 
?7max oc m^- 5 for v = 2. 

We have shown in this Section that we can use the com- 
bined constraints on the WIMP relic density and antipro- 
ton production to put an upper limit on r\, i.e. we can 
constrain the Hubble expansion at a very early time in 
the history of the Universe, as rj = A(Tf). The exact tem- 
perature at which we constrain the expansion in our dis- 
cussion, i.e. the standard freeze out temperature Tf, does 
however depend on the WIMP mass. Approximately, X{ 
is always of the order of 25 (xf = m x /T{). Some nu- 
merical examples are: Tf(m x = 60GcV) ~ 2.6 GeV, 
T f (m x = 200 GeV) ~ 7.5 GeV, T t (m x = 600 GeV) ~ 20 
GeV and T f (m x = 1000 GeV) ~ 32 GeV, where we have 
used the upper antiproton bound for the value of the 
WIMP annihilation cross section. 

Fig. [7| shows the upper bound on 77 as a function of 
T re for two different values of the WIMP mass and four 
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different values of the parameter v. Models that stay 
below the upper bound of 77 would give WIMPs that 
either are under-abundant or produce less antiprotons 
than the maximally allowed flux. Fig. [7| shows that the 
upper bound on r\ is independent of T ro as long as the 
latter is much smaller than the freeze out temperature. 
When the "re-entering" temperature is around one or 
two order of magnitudes below the standard freeze out 
temperature, ?7 max starts to increase and then becomes 
practically unbounded from above. In other words, one 
can increase the Hubble rate arbitrarily if this is done 
for a short amount of time. We see from fig. 0that the 
"re-entering" temperature at which 77 starts to get un- 
bounded is lower for the higher mass case despite the fact 
that the freeze out temperature is higher for m x = 1000 
GeV than for m x = 200 GeV. This can be understood 
from the analytic solution of the modified Boltzmann 
equation. We compare the contribution from the integra- 
tion from the modified freeze out until the "re-entering" 
temperature (i.e. the time where we still have enhanced 
expansion) with the contribution from the integration 
from the "re-entering" temperature untill today. The 
first contribution dominates when the "re-entering" tem- 
perature is very low. The analytic analysis shows that 
the two contributions becomes equal approximately when 
(Tf-T re )/T Ie = r)[h cS (x IC )yJg cS (xi)]/[h oS (xf) y/ g cS (x IC )} , 
where Tf is the freeze out temperature in the modified 
scenario and where we have set v = for definiteness. If 
77 stays constant then flh 2 — > (f2/i 2 ) st d when (Tf — T rc ) 
shrinks to zero. Alternatively, to avoid the WIMP to 
become under-abundant we have to increase 77 indefi- 
nitely. Since T grows approximately linearly with mass 
(Tf = m x /x{ and Xf depends logarithmically on the mass 
and the cross section) while the upper bound on 77 (fig. 
|SJ) is approximately proportional to the mass squared, 
it follows that r\ becomes unbounded at lower T rc for 
m x = 1000 GeV compared to m x = 200 GeV. 

The constraint on 77 = A(T[) immediately constrains 
the Hubble enhancement function A(T) given in Eq. (J5J 
at any temperature once we know the WIMP mass and 
the cosmological model (i.e. the parameters v and T re ). 
In Figs. 00 El E]| and 1111 we show the maximum allowed 
enhancement function for different values of the WIMP 
mass and the v parameter. The curves are derived for 
T re = 10~ 3 GeV. As we just saw, the upper bound on 77 
does not depend on T re when T IC <C Tf. Thus, a change 
of the re-entering temperature would solely change the 
temperature at which A(T) -» 1. Figure 151 151 and HP 
may be directly applied to constrain cosmological models 
with enhanced Hubble rate. For instance, for a scalar- 
tensor model like the one studied in Ref . [4J , which refers 
to a situation close to y = —1 (the actual value was 
v = -0.82) and T re = 0.1 GeV, by comparing Fig. EJhere 
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FIG. 8: Maximum allowed values of the enhancement function 
A(T) as a function of the temperature for our parametrization 
with v = — 1 (close to a scalar-tensor cosmology). Notice that 
time is running from right to left. The limit is derived from 
77max- The solid red line refers to m x = 100 GeV, the dashed 
green line to m x = 400 GeV, the dash-dotted blue line to 
m x = 600 GeV and the dotted purple line to m x = 1 TeV. 
For all the curves we use T le = 10 -3 GeV. The standard 
WIMP freeze-out temperature is given by the WIMP mass 
and the antiproton bound for the WIMP annihilation cross 
section. 

and Fig. 6 of Ref. we can conclude that the model 
of Ref. 0] is compatible with antiproton data only if the 
dark matter is composed by a WIMP with mass larger 
than 1 TeV. This example shows how constraining the 
antiproton data may be on otherwise viable dark energy 
models. 



B. The enhancement of the relic density 

In this section we derive upper limits on the enhance- 
ment of the WIMP relic abundance in modified cosmolo- 
gies as compared to the standard case, under the con- 
straint that fi/i 2 in the modified scenario satisfies the 
bound of Eq. |JTJ . This constraint alone is not actually 
enough to constrain the enhancement of the density, as 
we shall see. When the antiproton bound for the WIMPs 
is applied we find an upper bound on the density enhance- 
ment which is practically independent of the parameter 
v, but that is of the same order of magnitude as the upper 
bound on 77. 

As in the previous section we start by making a study 
in the (77, v) plane for a WIMP mass m x = 500 GeV. As 
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FIG. 9: Same as Figure[S]but for v = (overall enhancement 
of Hubble rate, with standard temperature evolution). 
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FIG. 11: Same as figure but for v = 2 (RSII model UU of 
extra-dimension cosmology) 
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FIG. 10: Same as figure|H]but for v = 1 (kination case). 

for Fig. 0] we make a scan in (cr ann w) for each (77, f) and 
keep models where the WIMP density, as calculated with 
the enhanced Hubble rate, fulfills the density constraint. 
For each of these models we calculate also the WIMP 
density as it would be in the standard case and we then 
find the ratio R = (nh 2 )/(nh 2 ) GR . 

Fig. 1 121 shows contours of values of R. Going from low 
to high values of r/, the different regions correspond to 
1 < R < 10, 10 < R < 100, 100 < R < 1000, 1000 < R 
with the highest value being around 7.5 • 10 3 . If we had 
continued to higher values of rj we would have obtained 



even larger enhancements of the WIMP density. From 
the previous section we know that rj can be bounded from 
above by requiring that the WIMP fulfills the constraint 
coming from the cosmic antiproton data. For m x = 500 
GeV the bound on 77 was found in Fig. 0] When we 
compare this figure with the contours of R in fig. 1121 
we conclude that the antiproton data leads to the upper 
limit R ~ 100 for m x = 500 GeV. We also see that 
the shape of the contours of R are very similar to the 
band in Fig. 0] which marks the limiting cross section. 
We therefore expect the upper bound on R to be almost 
independent of v. The final thing to note is that the 
numerical value of the upper bound on R and 77 are of 
the same order of magnitude. 

The upper bound on the relic density enhancement R 
as a function of the WIMP mass is shown in fig. ^3 for 
the four representative values of the parameter v. The 
approximate behaviour R <x rn^ 9 can be understood by 
means of the analytical approximation we already dis- 
cussed. 



C. Changing the assumptions 

In this Section we show how the upper bound on 
the Hubble enhancement 77 moves when some of our as- 
sumptions are changed. We are going to consider three 
kinds of modifications: we shall consider a temperature- 
dependent annihilation thermal cross section, by adding 
a p-wave contribution; we change the parametrization of 
the Hubble enhancement A(T); finally we take into ac- 
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FIG. 12: Contour plot of the enhancement R — 
(flh 2 )/ (Q,h 2 )Gn of the WIMP relic abundance in a scenario 
with enhanced Hubble rate compared to the standard GR 
cosmology. The different bands refer to (from left to right): 
1 < R < 10, 10 < R < 100, 100 < 7? < 1000, 1000 < R. The 
highest value of R is around 7.5 • 10 3 . We have fixed m x — 500 
GeV and T rc = 10" 3 GeV. For all points, the WIMP relic- 
density, as calculated in the modified cosmology, satisfies the 
dark matter density constraint. 



count the uncertainty of the propagation and diffusion of 
the cosmic antiproton signal. 

Let us first consider what happens if our assump- 
tion about pure s-wave annihilation does not hold. In 



this case (<r a 



relevant for the relic abundance, and 



(cann^)o, relevant for the antiproton calculation, are no 
longer equal. By adding a p-wave term b/x in the tem- 
perature expansion of Eq. (|19fl . (er ann u) turns out to be 
larger (or smaller, depending on the sign of b) in the early 
Universe than it is today where the p-wave contribution 
is negligible. The ensuing limit we obtain is less (more) 
stringent than the one we would obtain for pure s-wave 
annihilation, depending whether b is positive or nega- 
tive. Fig.^Jshows the effect for two representative cases 
in terms of r p = b/(ax{): r p = 0.5 (lower dash-dotted 
purple line) and r p = 10 (upper dash-dotted green line). 
The last value is actually extreme and it has been con- 
sidered to make an example. 

We notice that the results we obtain for a non negligi- 
ble p-wave case are similar to what we would obtain by 
assuming a value for the ratio 7Z in Eq. I|20|l of the order 
of TZ = (1 + u> r p ), where u> ~ 0(1) depends on details of 
the integrals given in Eq. I|16|) , and therefore also on the 
cosmological model (simple analytical expressions may 
be derived by means of the approximations of Appendix 
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FIG. 13: Upper limit of R = (Qh 2 ) / (f2/i 2 ) G R as a function of 
the WIMP mass. The bound comes from the combination of 
the WIMP relic density constraint and the antiproton bound. 
The lines, which practically fall on top of each-other, refer to: 
v = — 1, v = 0, v = 1 and v — 2. The mass lower limit de- 
pend on v (and it depends on the assumption of astrophysical 
parameters, which have been fixed to their median set in the 
figure). In the figure we have fixed T rc = 10 -3 GeV. 



|A"|) . The cases shown in Fig. I14lrefer to 7Z ~ 1.2 — 1.5 for 

r p = 0.5 and 7Z ~ 6— 10 for r p = 10. The figure therefore 
shows that the bound on 77 max is sensitive to large values 
of the ratio 71, as expected. In the case, for instance, of 
coannihilations, for which 1Z may be a large number (even 
some orders of magnitude) we see that the bounds from 
antiprotons are much less stringent. Approximately, the 
change in the upper bound on rj obtained in Eqs. for 
the case of time-independent cross section, can now be 
rephrased as: 



TL (o-annt')o (^ 2 )c7.DM 



'In 



C(y — 1) Xf 

TZ (o- a nn^)o (^h 2 ) C DM 

C(l — v ) X[ 



for v > 1 
for v < 1 



with, again, a slightly more complicated expression for 
u=l. 

Now, let us go back to the situation where we have only 
a temperature-independent thermal cross section but in- 
stead we modify the enhancement of the Hubble rate. 
Let us assume, instead of the enhancement A(T) defined 
in Eq. (J5J, the following function: 



A 2 {T) = 1 + jj 



T 

— ) tanh 2 



T-T„ 
aT I0 



(26) 
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FIG. 14: Upper bound on 7/ as a function of the WIMP mass 
under different assumptions, as derived from the combined 
constraint on the WIMP relic abundance and antiproton flux. 
All curves have T rc = 10 -3 GeV and v = — 1. The reference 
case is the red solid curve which is identical to our previous 
result in Fig. The uppermost and lowermost dotted lines 
delimit the band due to the uncertainty on the astrophysi- 
cal parameters for the antiproton propagation and diffusion. 
The uppermost curve is for the minimal set of astrophysi- 
cal parameters while the lowermost is for the maximal set of 
parameters. The two dash-dotted lines, instead, include a p- 
wave contribution for the annihilation in the early Universe. 
For the upper green curve the p-wave factor is r p — 10 (see 
text) and for the lower purple curve we have r p — 0.5. The 
dotted line which falls on top of the red solid (mean) line is 
calculated with the Hubble enhancement function A2 (T) with 
a change in the post-freeze-out temperature behaviour given 
by the parameter a = 100. 

The difference between A(T) and A^iT) is that in the 
latter case the tank is squared and in its argument we 
have introduced a factor a. This modification aims at 
obtaining (for large values of a) a "re-entering" into the 
standard regime less steep than in the previous case. This 
is illustrated in fig. ITS! 

A comparison between the upper bound on r\ in the 
case of a Hubble enhancement A(T) and A2(T) is shown 
again in Fig.^] The central blue dotted line is for A2 (T) 
with a = 100, T rc = 10 -3 and v = — 1. It can be seen 
that this line practically falls on top of the limit (red solid 
line) derived for A(T) (and T rc = 10~ 3 , v = -1). This 
result shows that the actual shape of the enhancement 
function around "re-entering" is not of major importance 
as long as the "re-entering" occurs long after the WIMP 
freeze out. 

Let us finally consider how the upper bound on r\ is 
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FIG. 15: Enhancement function ^(T 1 ) of Eq. 126H as a func- 
tion of the temperature. Notice that time runs from right to 
left. The solid red line has v = —1, the dashed green line 
v = 0, the dash-dotted blue v = 1 and the dotted purple line 
has v — 2. We have chosen a — 100. All other parameters are 
identical to those of A(T) in Fig. □ i.e. 77 = 10 2 , T ro = 10 -2 
GeV and T t = 17.3 GeV. 

affected by the uncertainty in the astrophysical param- 
eters of the propagation and diffusion of the antiproton 
flux. We refer for this uncertainty to Ref. |l2|. Up 
to now, we have used a set of astrophysical parameters 
which provides the median value of the antiproton flux 
[l3 | . This set of parameters is the one which best fits cos- 
mic ray observations, mainly the B/C ratio j2||. We now 
show in Fig.^lth results obtained for the allowed range 
of astrophysical parameters, which give a maximal and 
minimal antiproton flux. The uppermost dotted curve in 
Fig. ^] shows the result obtained using the "minimal" 
set of astrophysical parameters for the propagation and 
diffusion, while the lowermost dotted curve shows the re- 
sult as calculated using the "maximal" set. These curves 
should be compared to the central solid curve, which is 
our previous result obtained with the best-fit parame- 
ters. For all the curves we have T re = 10 -3 GeV and 
v = — 1. We see that the uncertainty in the propaga- 
tion function causes an uncertainty of almost one order 
of magnitude in both direction for the upper limit of r), 
as expected from the discussion in Sec. [V]and from the 
analytical approximations of Eq. (|25ll . 

From Fig. ^]we also see that a different choice on the 
astrophysical parameters may lead to a different lower 
limit on the WIMP mass, as a consequence of the fact 
that lowering the bound on the cross section (for the 
maximal set of astrophysical parameters) causes the up- 
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per bound on 77 to decrease. This can be seen from the 
plot of the (r), v) plane in Fig. EI For a selection of masses 
this shows the bound as calculated with the median prop- 
agation function. For a given mass, the region to the left 
(i.e. at smaller 77) is allowed by the antiproton data as the 
cross section here is lower. Thus, when we decrease the 
upper bound on the cross section, i.e. increase the effect 
of the propagation parameters, all the bounds in Fig. 
will move to smaller 7/ (in particular the lower masses 
will move outside the plot). This is consistent with the 
results of Ref. El. 



VII. ANTIPROTON FLUX CONSTRAINTS ON 
SOME SPECIFIC MODELS 



We know from Sec. I VI 61 that 77 is approximately equal 
to the WIMP density enhancement. This is in agreement 
with the conclusions of the earlier references on kination 
and enhanced WIMP density; see e.g. Ref. [f| and Ref. 
0. In Ref. they found that the maximal density 
enhancement compatible with the BBN bounds is of the 
order of 10 6 (for WIMP masses smaller than 1 TeV), i.e. 
?/$ < 10 12 . With our argument based on the maximum 
allowed antiproton flux we can derive a much stronger 
constraint on this scenario. The most conservative bound 
is found for large masses. For m x < ITeV (and v = 
1) we have 77 < 0(1O 3 ) and so 77$ < 10 6 . While for 
light WIMPs (m x < O(10 2 GeV)) we have the much more 
stringent bound 77 < 10 and so 77$ < 10 2 . 



As we discussed in the previous sections, the maxi- 
mum allowed antiproton flux gives interesting model- 
independent constraints on the maximum allowed en- 
hancement of the expansion rate at dark matter freeze 
out. Let us now discuss these constraints in connection 
to the specific cosmological models we shortly introduced 
in Section liTTl 

A. Scalar-tensor theories 



As mentioned in section llll Al we do not have a general 
analytical correspondence between the Hubble expansion 
rate in the scalar-tensor model and our parametrization 
(JSJ). However, the numeric example of reference which 
assumes a slowly varying scalar field, can be described by 
our parametrization with v = —0.82 and 77 = A(ip(Ti)), 
where A(<p) is the Weyl factor, which in this case gives 
the Hubble enhancement function A(T). From figure 
we get the upper bound: 

A(tp(T { )) < O(10 2 ) for m x < O(500GeV) (27) 

and Tf = O(lOGcV). From figure |H1 we can also find 
the upper bound on A(T) = A(ip) at earlier times. Wc 
see that A(ip) < 10 for T > O(10 2 GeV) and m x < 
O(500GeV). It is interesting to observe that our re- 
sults strongly constrain Fig. 5 of in which it could 
seem that properly choosing the initial conditions for the 
scalar field, each enhancement for the expansion rate is 
allowed. Antiproton data limit this possibility. 



C. Extra Dimensions 

As mentioned in section IIII CI the RSII model [ijj 
gives rise to an enhanced expansion rate which can be 
described by our parametrization with the values of the 
parameters v = 2 and 77 = */ p r (Tf) / (2A). Where A is 
the tension on the brane and p r is the radiation ener- 
gy density. The upper bound on 77 therefore gives us 
a lower bound on the brane tension. The upper bound 
on 77 for v = 2 is shown in Fig. [fJl For 77 < O(10 2 ), 
which is true for WIMP masses m x < £>(10 2 GeV), we 
get the most stringent bound. With Tf ~ 5 GeV and 
g cS (T f ) ~ 9 2 wc find the bound A > 0(1 GeV 4 ). A 
more conservative bound is found for higher masses. For 
m x = 0(500 GeV) we have 77 < O(10 4 ) and T t ~ 20 GeV, 
leading to A > O(10- 2 GeV 4 ). Finally using the relation 

El 



A = 



M| 



4ttM 2 



(28) 



where M p i is the four dimensional Plank mass and M5 
the five dimensional one, we can derive a lower bound 
on M5. For the most stringent bound on A we get 
M5 > 3 • 10 3 TeV, while the more conservative bound just 
give a small change, M5 > 10 3 TeV. We observe that the 
BBN bound is A > 1 McV 4 and so M 5 > 30 TeV, 0. 
Therefore in this context, our analysis based on the an- 
tiproton flux is much more stringent than the one coming 
form the BBN constraint. 



B. Kination 



VIII. CONCLUSION 



In sec. IIII Bl wc found that v = 1 and 77 = y/rj$ = 
y/ p$ (Tf ) / p r (Tf ) can be used to describe the enhanced 
expansion rate during a kination phase. 



The expansion rate of the Universe in the era before 
BBN is predicted to be different from the one of standard 
FRW cosmology in a host of cosmological models which 
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are derived either from modification of General Relativ- 
ity, or from brane physics, or from attempts to provide a 
consistent explanation to dark energy. In the present pa- 
per we concentrated on models which predict an enhance- 
ment of the expansion rate, like e.g. those discussed in 
Ref. where the dark energy problem finds a solution 
in a scalar-tensor theory of gravity with an exponential 
run-away behaviour of a scalar field, in kination models 
like those of Ref. 0, EH Rj Q , in D-brane models like the 
RSII model of Ref. |llj or models of modified cosmol- 
ogy like those in Ref. 0, . These enhancements may 
be sizeable, without evading the post-BBN bounds (light 
primordial element production, CMB anisotropics, grav- 
itational tests) which may be applied to these cosmologi- 
cal models: the enhancements may reach up to four 0] or 
even six 0,0 orders of magnitude. However, an imprint 
of this enhancement may be left on the properties of the 
particles responsible of the dark matter. This happens 
if the enhancement of the expansion rate occurs close to 
the time when the dark matter particles decouples from 
the thermal bath: this possibility relies on the fact that 
the cosmological bound on the amount of dark matter 
can be fulfilled by particles with larger annihilation cross 
section, as compared to standard cosmology 0, |j| , due 
to an anticipated freeze-out. In this situations, indirect 
detection signals are typically enhanced, as long as the 
thermal cross section responsible of the relic abundance is 
close to the one which determines the indirect detection 
signals (like e.g. when the typical temperature expansion 
of the thermal cross section is applicable). 

Cosmic antiprotons produced by WIMPs in the Galaxy 
are currently the best option for constraining particle 
dark matter annihilation cross sections 0, Il3l| . since 
they are not strongly dependent on assumptions on the 
dark matter density profile, like instead is the case of the 
gamma-ray signal. Also, cosmic antiprotons produced by 
WIMPs are naturally predicted to be at the level of the 
background and of the data, when the WIMP accounts 
for the CDM content of the Universe, and without invok- 
ing high degrees of over-densities in the galactic neigh- 
bourhood [13I . Since the experimental data are in good 
agreement with the expected background (3, I22}, not 
much room is left for an exotic component and therefore 
antiprotons may be used to set constraints [l^.lT^|. 

In this paper we exploited the effect induced on the 
antiproton signal by an increased expansion rate, to set 
bounds on the maximal enhancement of the expansion 
rate which can be allowed in the early Universe in a pre- 
BBN epoch Q. These limits apply to an evolutionary 
phase close to dark matter decoupling (i.e. at tempera- 
tures of the Universe in the range of 1-30 GeV, for dark 
matter of about 10-1000 GeV mass) and are subject to 
the (quite reasonable) hypothesis that dark matter is a 



thermal relic. The bounds we have obtained are quite 
strong, especially for light dark matter. Our main result 
is summarized by Fig. El F° r WIMP masses lighter than 
100 GeV, the bound on the Hubble-rate enhancement 
ranges from a factor of a few to a factor of 30, depending 
on the actual cosmological model, while for a mass of 500 
GeV the bound falls in the range 50-500. These bounds 
are affected by the still large astrophysical uncertainties 
in the calculation of the antiproton signal |l2j , which re- 
flects in a factor of 5 more stringent, or a factor of 10 
more loose bounds. Nevertheless, the limits we obtain 
are much more stringent than the possible enhancement 
of these cosmological models. 

A caveat is in order here: whenever the effective ther- 
mal annihilation cross section in the early Universe (re- 
sponsible for the relic abundance) is sizably larger than 
the one in the galactic halo (which determines the an- 
tiproton signal) the bounds we determine are necessarily 
loosened. This occurs when coannihilation effects are 
relevant in determining the relic abundance. In this case 
the annihilation cross section in the Galaxy is typically 
much smaller than the effective one which sets the de- 
coupling: in this case our limits still apply, but they are 
less constraining, approximately by the ratio of the two 
cross sections. We nevertheless notice that the occur- 
rence of coannihilation is usually an accidental feature. 
Another case where our limits are less severe is when 
WIMP annihilation occurs mostly into leptons: in this 
case antiprotons are not produced and our bounds are 
loosened by a factor given by the branching ratio into 
non-lcptonic final states. In all other situations, which 
represent the most typical cases, the limits we quote are 
actually necessarily present. 



Finally, we applied our bounds to set constraints on 
some specific models. In the case of the scalar-tensor 
dark energy models of Ref. 0, we constrain the en- 
hancement to be less than a factor 100 for dark matter 
lighter than 500 GeV. In the case of kination models, the 
maximal enhancement is 1000 for WIMPs lighter than 
1 TeV, while for dark matter lighter than 100 GeV the 
enhancement must be below 10. Finally, for the Randall- 
Sundrum D-brane model of Ref. |ll| . our limit may 
be used to set bounds on the tension on the brane: for 
WIMPs lighter than 100 GeV, the string tension must be 
larger than 1 GeV 4 ; for heavier dark matter, the bound 
is 10 ~ 2 GeV 4 . These limits reflect in a lower bound of 
about 10 3 TeV on the 5-dimensional Planck mass and 
are much more severe than the usual constraints derived 
from BBN physics. 
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APPENDIX A: APPROXIMATE 
BOLTZMANN-EQUATION SOLUTIONS FOR 
THE MODIFIED HUBBLE RATE 

The Boltzmann equation Eq. (|12|) may be easily solved 
for a modified Hubble rate as defined in Eq. J5J) by suit- 
able approximations. We first make explicit the tempera- 
ture (or x = ni x /T) dependence in Eq. 1|12|) by rewriting 
the equation in the usual form: 

dY 
dx 



h ef f(x) 



-(a ann v)(Y 2 -Y 2 ) 



45G x /A^(x)g cS (x) x< 

(Al) 

where G is the Newton constant. The standard cosmo- 
logy case is recovered for A(T) = 1. 

We first define the temperature of particle freeze-out 
as the temperature Tf when: 



Y{x { ) = (1 + c)y eq (z f ) 



(A2) 



with c a constant of order 1 and Y cq (x) ~ 45 x 
2~ 5 / 2 7r~ 7 / 2 g h~g (x) x~ 3 / 2 exp(— a;) is the equilibrium 
abundance of a non-relativistic particle with internal de- 
gree of freedom g. Condition (|A2() and Eq. IjAljl lead to 
the following implicit expression for the freeze-out tem- 
perature: 



Xf 



In 



2c 



1 



-0.038 rripigm^ 



x) /2 A{xf)gl£(xf) 



(A3) 



where m p \ is the Planck mass. For definitencss, we use 
c = 1 throughout the paper. The standard freeze-out 
temperature is obtained when A(T) = 1 around the time 
of particle decoupling. The models we are considering, 
where A(T) > 1 at that epoch, imply a smaller Xf, i.e. a 
larger freeze-out temperature Tf. 

The current value Yq of the comoving abundance is 
then obtained by integrating Eq. (|A1() from Xf down to 
xq — > co, neglecting the WIMP production term: 



dY 
dx 



h eS (x) 



■(a wn v)Y 2 (A4) 



■ dx 



(A5) 



45G y/A*(x)g og (x) x 2 
which gives the solution: 

Y Y cq (x f ) x J xt A(x): 

where B = (n/A5G)^ 2 and Q(x) = h cS (x) / y/ g cS (x) . 
The WIMP relic abundance is then obtained by means 
ofEq. 

For large enhancements of the Hubble rate around the 
WIMP freeze-out, i.e. for large values of the parameter 
•q, we may approximate A{T) as a power-law function 
with exponent v. If v < 0, we may assume this behaviour 
all the way down to the re-entering temperature x ra - 



a(t) ~ n 



(rjx^)x' 



(A6) 



When v > the above approximation may not be used 
up to x re , because A{x) approaches unity as the Universe 
cools down (see Fig. ^) and therefore we cannot neglect 
the 1 in Eq. (JSJ). We define as X\ the temperature when 

(T) Xf) x~ 



1, i.e. x\ = XfT] 1 ^. In this case we approx- 



imate: 



A(T) 



(r)x\)x- v 
1 



X < X\ 



(A7) 



Below a; rc , in any case A(T) = 1 by definition. 

With these approximations, the solution of Eq. (|A5() 
for a temperature independent annihilation cross section 



a is: 



1 



1 

% ~ Y cq ( Xi ) 

+ Bm x aQ(xf) 



(A8) 



T(x{,x,rj, v) + r G - 

x 



where x = min[xi,x rc ] and tq = Q{x)jQ(xi) is typically 
around 0.5-1 for WIMPs masses in the range GeV-TeV. 
For large enhancements (r\ 3> 1) and a low re-entering 
temperature, typically: x = x\ for v > and x = x re for 
v < 0. In Eq. (|A8|) we have assumed a slowly-varying 
Q{x) in the integrals. 

In our discussion in the text, typically the relevant 
term in Eq. (|A8|) is the one which contains the func- 
tion X(xi , x, ?y, v): the other two terms may be usually 
neglected to a few percent level of approximation. This 
relevant function is easily found to be: 



1 



1 



v — 1 Xf r\ 1 l v 



for v > 1, then 



1 = 



1 



V Xf J] 



for v < 1 and finally 



I 



1 Info) 
X[ 77 



(A9) 



(A10) 



(All) 



for v = 1 (kination). 



ACKNOWLEDGMENTS 

We acknowledge Research Grants funded jointly by the 
Italian Ministero dell'Istruzione, dell'Universita e della 
Ricerca (MIUR), by the University of Torino and by the 
Istituto Nazionale di Fisica Nucleare (INFN) within the 
Astroparticle Physics Project. R. Catena acknowledges 
a Research Grant funded by the VIPAC Institute. 



19 



[1] D.N. Spergel et al . (WMAP Collaboration), 
|arXiv:astro- ph/0603449] (submitted to Astrophys. 
J.) 

[2] A reduction of the expansion rate is also possible; see 
e.g. Ref. Q. We are not going to discuss these models 
here, since when the Hubble rate is rduced, the antipro- 
ton bounds derived in the paper does not lead to relevant 
constraints. 

[3] R. Catena, M. Pietroni and L. Scarabello, Phys. Rev. D 

70, 103526 (2004) [arXiv:astro-ph /0407646] ; 

C. de Rham, T. Shiromizu and A. J. Tolley, 

arXiv:gr-qc/0604071 . 
[4] R. Catena, N. Fornengo, A. Masiero, M. Pietroni 

and F. Rosati, Phys. Rev. D 70, 063519 (2004) 

|arXiv:astro-p h/04036i4|. 
[5] P. Salati, Phys. Lett. B 571, 121 (2003). 
[6] F. Rosati, P hys. Lett. B 570, 5 (2003) 

arXiv:hep-ph/0309124 . 
[7] S. Profumo and P. U llio, JCAP 0311, 006 (2003) 

|arXiv:hep-ph/0309220] 
[8] C. Pallis, JCA P 0510, 015 (2005) 

|arXiv:hep-ph/05030'80] . 
[9] J. D. Barrow, Nucl. Phys. B 208, 501 (1982). 
[10] M. Kamionkowski and M. S. Turner, Phys. Rev. D 42, 

3310 (1990). 

[11] L. Randall and R. Sundrum, Phys. Rev. Lett. 83, 4690 
(1999). 

[12] F. Donato, N. Fornengo, D. Maurin, P. Salati 

and R. Taillet, Phys. Rev. D69. 063501 (2004) 

|arXiv:astao~p h/03062d7]. 
[13] A. Bottino, F. Donato. N. Fornengo and P. Salati, Phys. 

Rev. D72, 083518 (2005). 
[14] R. H. Cyburt, B. D. Fields, K. A. Olive and 

E. Skillman, Astropa rt. Phys. 23, 313 (2005) 

|arXiv:astro-ph/0408033"l . 

[15] Y. Wang and P. Mukherjee, [arXiv:astr o-ph /060 405 1 1 . 
[16] G. B. Zhao, J. Q. X ia, B. Feng and X. Zhang, 

|arXiv:astro-p h/060362l"1 . 
[17] J. Albert et al. [SNAP Collaboration], 

|arXiv:astro-ph/0507460"l . 
[18] B. Bertotti, L. less and P. Tortora, Nature 425, 374 

(2003). 

[19] T. Damour, arXiv:gr-qc/9606079 , lectures given at Les 

Houches 1992, SUSY05 and Corfu' 1995. 
[20] G. Esposito-Farese and D. Polarski, Phys. Rev. D 63, 

063504 (2001) [arXiv:gr -qc /0009034| . 
[21] R. H. Dicke, Phys. Rev. 125, 2163 (1962). 
[22] R. Catena, M. Pietroni and L. Scarabello, 

|arXiv:astro-ph/0604492| . 
[23] P. J. Steinhardt, L. M. Wang and L Zlatev, Phys. Rev. D 

59, 123504 (1999) |arXiv:astro-ph/9812313| . 
[24] T. Shiromizu, K. I. Maeda and M. Sasaki, Phys. Rev. D 



62, 024012 (2000). 

[25] R. Durrer, AIP C onf. Proc. 782, 202 (2005) 
[arXiv:hep-th/0507006l . 

[26] F. Donato, D. Maurin, P. Salati, A. Barrau, G. Boudoul, 
and R. Taillet, Astrophys. J. 563, 172 (2001). 

[27] D. Maurin, R. Taillet, F. Donato, P. Salati, A. Barrau, 
and G. Boudoul, Research Signposts (2004), "Recent Re- 
search Developments in Astronomy and Astrophysics", 
v.2, p. 193 [arXiv:astro-ph /02121lT| . 

[28] D. Maurin, F. Donato, R. Taillet, and P. Salati, Astro- 
phys. J. 555, 585 (2001). 

[29] S. Orito, et al. (BESS Collaboration), Phys. Rev. Lett. 
84, 1078 (2000). 

[30] T. Maeno, et al. (BESS Collaboration), Astropart. Phys. 
16, 121 (2001). 

[31] M Aguilar, et al. (AMS Collaboration), Phys. Rep. 366, 
331 (2002). 

[32] M. Boezio, et al. (CAPRICE Collaboration), Astrophys. 

J. 561, 787 (2001). 
[33] J. Lavalle, J. Pochon, P. Salati and R. Taillet, 

[arXiv:astro-ph /0603796| . 
[34] F. Donato, N. Fornengo, P. Salati, Phys. Rev. D62, 

043003 (2000). 

[35] K. Mori et al, Astr ophys. J. 566, 604 (2002) 

[arXiv:astro-p h /0109463"1 . 
[36] M. Aguilar et al. (AMS Collaboration), Phys. Rep. 

366/6, 3 31 (2002). 

[37] P. Ullio, [arXiv:astro- ph /9904086] . 

[38] S. Mizuta and M. Yamaguchi, Phys. Lett. B 298, 
120 (1993); J. Edsjo and P. Gondolo, Phys. Rev. D 
56, 1879 (1997) [arXiv:hep-ph/9704361| ; J. R. Ellis, 
T. Falk and K. A. Oliv e, Phys. Lett. B 444, 367 (1998) 
|arXiv:hep-ph/9810360] . 

[39] P. Binetruy, G. Girardi and P. Salati, Nucl. Phys. B 237. 
285 (1984) ; K. Griest and D. Seckel, Phys. Rev. D 43, 
3191 (1991) 

[40] P. Gondolo, and G. Gelmini, Nucl. Phys. B 360, 145 
(1991). 

[41] This is the case when the cross section does not exhibit 
a strong temperature-dependence and when thermal ef- 
fects like coannihilations are not effective in the determi- 
nation of the thermal average of the cross section. This 
is the case we shall focus on mainly in this paper. In 
Scc lVI 01 we study a temperature dependent cross sec- 
tion. 

[42] Note that models that have the same WIMP mass and 
cross section would have identical relic abundance in 
standard cosmology. Here, however, we are interested in 
the relic abundance in a cosmology with enhanced Hub- 
ble rate, and this abundance does depend also on the 
parameters r/ and v. 



